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Abstract 



A search for charged Higgs bosons produced in pairs is performed with data col- 
lected at centre-of-mass energies ranging from 189 to 209 GeV by ALEPH at LEP, cor- 
responding to a total luminosity of 629 pb -1 . The three final states t + v t t~u t , cst~v t 
and esse are considered. No evidence for a signal is found and lower limits are set on 
the mass mjj± as a function of the branching fraction B(H + -^t + v t ). In the framework 
of a two-Higgs-doublet model, and assuming B(H + -^t + ^ t )+B(H + — ►cs)=l, charged 
Higgs bosons with masses below 79.3GeV/c 2 are excluded at 95% confidence level 
independently of the branching ratios. 
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1 Introduction 



The Standard Model of electroweak interactions requires only one doublet of complex scalar 
fields, resulting in a single neutral Higgs particle. The simplest extensions of the Standard 
Model assume two complex scalar-field doublets, with a total of eight degrees of freedom. As 
in the Standard Model, three of the degrees of freedom are associated with the longitudinal 
components of the W and Z bosons. The remaining five degrees of freedom appear as five 
physical scalar Higgs states: three neutral Higgs bosons and the charged Higgs bosons H 1 * 1 . 

In the two-Higgs-doublet case, the charged Higgs boson couplings are completely 
specified in terms of the electric charge and the weak mixing angle #w The production 
cross-section thus depends only on the mass mn± • For masses accessible at LEP 2 energies, 
the charged Higgs boson decays with negligible lifetime and width into either cs/cb or r + u T 
final states. Because the analyses are not sensitive to the quark flavour, and because the 
cs decay mode dominates over cb, cs stands for either cs or cb in the following. Therefore, 
B(H + ^t + z/ t )+B(H + ^cs)=1 is assumed and H + H pair production leads to three final 
states {j Jr v T T~v T , csr _ z/ T /csr + z/ r and esse) for which separate searches are performed. 

The ALEPH data collected at energies up to 189 GeV have already been analysed and 
the search results published in Refs. |], 0, |3|. The negative result of the search, under the 
hypotheses specified above, was translated into a lower limit on the H 1 * 1 mass of 65.5 GeV/c 2 
at 95% confidence level (C.L.). Results from other experiments are given in Ref. Pj. The 
present letter describes the search for pair-produced charged Higgs bosons using the data 
collected up to the end of data taking. An improved analysis has been designed for the fully 
leptonic channel. In the semileptonic search, the rejection of the W + W~ background has 
been refined with a method based on a combination of the charge-tagged boson production 
angle and a r polarization estimator. For the four-jet event selection, the linear discriminant 
analysis (LDA) has been re-optimized to account for the additional integrated luminosity 
collected at increased centre-of-mass energies. 

2 The ALEPH detector and event samples 

A complete and detailed description of the ALEPH detector and its performance, as well 
as of the standard reconstruction and analysis algorithms can be found in Refs. [^, ^j. Only 
those items relevant for the final states under study in this letter are summarized below. 

The trajectories of the charged particles (called charged tracks in the following) are 
measured with the central tracking system, formed by a silicon vertex detector, an inner 
drift chamber and a large time projection chamber, all immersed in the 1.5 T axial magnetic 
field from a superconducting solenoidal coil. Electrons and photons are identified in the 
electromagnetic calorimeter, a highly segmented sampling calorimeter placed between the 
tracking device and the coil. Muons are identified in the hadron calorimeter, a 1.2 m thick 
iron yoke instrumented with 23 layers of streamer tubes, surrounded with two double layers 
of muon chambers. Together with the luminometers, the hermetic calorimetric coverage 
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extends down to 34mrad of the beam axis. The missing energy and momentum from, 
e.g., tau charged Higgs boson decays, are determined with an energy- flow algorithm which 
combines particle identification, tracking and calorimetry information into a set of energy- 
flow particles, used in the present analyses. 

The data analysed in this letter were collected at LEP between 1998 and 2000 at e + e~ 
centre-of-mass energies ranging from 189 to 209 GeV, corresponding to a total integrated 
luminosity of 629pb _1 . The details for each sample are given in Table |I[ 

Table 1: Integrated luminosities, centre-of-mass energy ranges and mean centre-of-mass energy values for 
the data collected with the ALEPH detector from 1998 to 2000. 



Year 


Luminosity (pb L ) 


Energy range (GeV) 


(V5> (GeV) 


2000 


217.2 


204 - 209 


206.1 


1999 


42.0 




201.6 




86.3 




199.5 




79.8 




195.5 




28.9 




191.6 


1998 


174.4 




188.6 



Fully simulated samples of events reconstructed with the same programs as the data 
were used for the background estimates, the design of the selections and the optimization 
of the selection cuts. The most important background sources are (i) difermion events 
(e + e _ — ► r + r~ and qq) simulated with the K0RALZ |7| generator; and (it) e + e~ — ► W + W~ 
and other four-fermion processes simulated with the K0RALW || and PYTHIA [[J generators. 
Event samples of these background processes, corresponding to at least 20 times the 
collected luminosity, were generated. The W + W~ cross sections predicted by RAC00NWW [[l(J 



and YFSWW |TTJ were used as discussed in Ref. il2| . Finally, the two-photon interactions 
(77 — *■ leptons) were simulated with the PH0T02 ]lj| generator. Samples of these events 
with at least six times the collected luminosity were generated. 

The signal events generated with the HZHA [|TJj] program were simulated for each of 
the final states and centre-of-mass energies (Table [I]), and for charged Higgs boson masses 
between 45 and 100GeV/c 2 . 



3 Analyses 

An event selection has been defined for each of the t + v t t~v t , cst~v t /cst + u t (hereafter 
referred to as csT~i> T ) and esse channels, and was optimized for B(H + ^t + u t ) = 100%, 
50% and 0%, respectively. The selection criteria were chosen to achieve the highest 95% 
confidence level expected limit on the charged Higgs boson mass in the absence of signal. 
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3.1 The t + v t t v T final state 



Events with two to six charged tracks (at least one and at most four of each sign) are 
considered. Leptonic events W + W~ — > tvt'v (£, £' = e or fi) are rejected by requiring that 
the momentum of any identified electron or muon be less than 0.1-y/i. The events are then 
forced to form two jets with the JADE algorithm [15]. An event is selected if both jet polar 
angles 9 12 satisfy | cos #1,2] < 0.96, if their reconstructed masses are less than 3 GeV/c 2 and 
if each jet contains at least one charged track. To suppress the high cross section 77 — > ff 
processes, the total visible mass is required to be in excess of 0.075-^/5, the momentum 
transverse to the beam is required to be greater than lOGeV/c, and there must be no 
energy deposited in a cone of 12° around the beam axis. The signal selection efficiency of 
the latter cut is corrected for the effect of the beam-related background, not included in the 
simulation, and is estimated from events triggered at random beam crossings. The relative 
loss of signal efficiency is about 7%. 

Nearly coplanar tau pairs from e + e~ — > t + t~{j) are rejected by requiring that the angle 
a between the two tau jets be less than 170° and the angle between the projections of their 
momenta onto the plane transverse to the beam axis be less than 165°. The missing energy 
is required to be greater than 80GeV and the missing mass greater than 70GeV/c 2 . In 
order to improve the W + W~ background rejection, an LDA has been used to construct a 
discriminant variable D from a combination of the following four quantities: 



• a charge-tagged angular variable calculated from the polar angles of the r jets and 
their charges as C = ~ [Qicos^i + <52Cos# 2 ]; 

• the angle a between the two tau jets; 

• the missing transverse momentum of the event P™ 1SS ; 

• the value 2/23 of the jet-clustering resolution parameter for which the transition from 
two to three jets occurs. 



The optimal discriminant variable was found to be 

D = 0.930 C - 0.250 a + 0.008 Pf ss - 110 y 23 + 0.426 , 

where a is in radians and P™ 1SS in GeV/c. The distribution of Dq is displayed in Fig. [1]. 
This quantity is used as a discriminant variable in the derivation of the mass limit. 

The signal event selection efficiencies, parametrized as a function of mn±, are given in 
Table for = 206 GeV. The selection efficiencies are almost independent of the centre- 
of-mass energy and increase only slightly with mn±. For a signal with mH±=85GeV/c 2 
and B(H + ^t + z/ t ) = 1, a total of 16.5 events is expected in the data taken at centre-of- 
mass energies between 189 GeV and 209 GeV. The numbers of events selected are given in 
Table [3], compared to the expectations from the Standard Model backgrounds, dominated 
by W + W~ production. 
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Figure 1: The distribution of the discriminant variable Dq described in the text for the fully- 
leptonic channel. The points are the data, the open histogram is the Standard Model background 
and the hatched histogram represents the Higgs signal expectation, absolutely normalized, with 
m H ± = 85GeV/c 2 . 



The systematic uncertainty on the number of expected signal events is estimated to be 
3.1%, dominated by the effect of limited Monte Carlo statistics (2.4%) and the uncertainty 
on the cross section for charged Higgs boson production (2%). The systematic error on the 
background level is estimated to be 1.5%, dominated by the effects of limited Monte Carlo 
statistics (1.3%), by the uncertainty on the cross section for the W + W~ process (0.5%) and 
the uncertainty on the cross section for two-photon production (5%). 



3.2 The csr v T final state 

The mixed final state cst~u t is characterized by two jets originating from the hadronic 
decay of one of the charged Higgs bosons and a r jet with missing energy due to the prompt 
neutrino as well as to the neutrino (s) from the subsequent r decay. 

The preselection is the same as that described in Ref. 0. In order to identify the r jet 



an algorithm based on "minijets" is used as described in Ref. pq . If a minijet satisfies the 
r-jet selection criteria, the rest of the event is clustered into two jets using the Durham [17] 
clustering algorithm. A kinematic fit is performed with the constraints of energy and 
momentum conservation and equality of the cs and t + u t masses. If there is more than one 
r candidate the combination with the lowest x 2 is taken. 
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Table 2: The signal event selection efficiencies e (in %), parametrized as a function of the charged 
Higgs boson mass m H ± , at y/s= 206 GeV. 



m H ±(GeV/c 2 ) 


60 


65 


70 


75 


80 


85 


90 


e (t + v t t~u t ) 


24.4 


25.5 


26.4 


27.3 


28.0 


28.5 


28.9 


e (csr~P T ) 


49.1 


48.0 


45.8 


42.8 


38.8 


33.9 


28.0 


e (esse) 


60.7 


62.9 


64.5 


65.5 


66.1 


66.3 


66.3 



Table 3: Numbers of candidate events and background expected from Standard Model processes, 
for each of the three years of data taking. 



Channel 




observed 


expected 




(GeV) 


events 


background 


T + V T T~U T 


188.6 


14 


11.0 




192-202 


22 


15.6 




204-209 


9 


14.0 


CST~P T 


188.6 


63 


67.3 




192-202 


89 


113.1 




204-209 


127 


108.9 


cssc 


188.6 


778 


826.3 




192-202 


1034 


1102.6 




204-209 


950 


963.2 



In order to reject background from W + W~ — > (e//i)z/qq', the measured energy of the 
r jet boosted into the Higgs rest frame is required to be less than 0.175a/s. The boost is 
performed using the information from the hadronic side of the event. 

After this procedure the following four variables are chosen to further suppress the 
background: 

• the total missing transverse momentum of the event, P™ ss ; 

• the isolation angle #i SO of the r, denned as the half-angle of the cone around the r jet 
direction containing 5% of the total energy of the rest of the event; 

• the x 2 from the kinematic fit; 

• the decay angle Of 1 , denned as the angle between the r momentum in the Higgs 
boson centre-of-mass frame and the Higgs boson flight direction, charge-tagged with 
the charge of the r, to exploit the asymmetry in the W system, absent for scalars. 

The four variables are linearly combined into one variable, D\, defined as 

D 1 = 0.021 Pf ss + 0.400 # iso - 0.058 X 2 ~ 0.148 9f - 0.881 
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Figure 2: (a) The distribution of the discriminant variable D 2 described in the text for the semi- 
leptonic channel, (b) The distribution of the fitted mass of the Higgs boson candidates after the 
cut on D 2 . The points are the data, the open histogram is the Standard Model background and 
the hatched histogram represents the Higgs boson signal expectation with m H ± = 75 GeV/c 2 . The 
signal is arbitrarily normalized. 



where P™ 1SS is in GeV/c, and 9- lso and 8 T h are in radians. Events are selected by requiring 
that D\ > —0.1. The background consists primarily of W + W~ — > £z/qq' events. 

Due to the scalar nature of the H + , the r + from its decay is produced in a left- 
handed helicity state, in contrast to the r +, s from W + decays. Variables designed for the 



measurement of the r polarization at LEP 1 [18] have been used to form an event-by-event 



helicity estimator, £ T . This variable, together with the charge-tagged production angle 
^prod 05 i s use d t° discriminate further between W + W~ — > ruqq' and H + H~ — > cst~T> t 
events. The two variables are combined into another variable, D 2 , defined as 

D 2 = - 0.461 0^ od - 0.517 S T + 1-020 , 

where 9^ od is expressed in radians. The distribution of D 2 is shown in Fig. 0a. The cut 
optimization yields D 2 > —0.3 for mn±=75 GeV/c 2 . The selection efficiencies are given in 
Table as a function of the Higgs boson mass for y/s = 206 GeV. They are only weakly 
dependent on yfs. In the data collected between y/s = 189 and 209 GeV, the numbers of 
selected events are compared with the background expectations in Table [| The fitted-mass 
distribution of the Higgs boson candidates is shown in Fig. [|b. For mH±=77GeV/c 2 , close 
to the sensitivity of this search, and for B(H + — >t + v t )=0.5, a total of 21.2 signal events is 
expected. 

The systematic uncertainty on the number of expected signal events is estimated to be 
3.0%. The main contributions are the finite size of the simulated event samples (2.2%), 
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calorimeter calibration uncertainties (0.5%) and the uncertainty on the cross section for 
charged Higgs boson production (2%). The systematic error on the background level was 
estimated to be 3.9%. The main contributions are from limited statistics of the simulated 
event samples (2.5%), uncertainty on the cross section for the W + W _ process (0.5%) and 
calibration uncertainties (3%). 



3.3 The esse final state 



The hadronic decays of pair-produced charged Higgs bosons lead to a four-jet final state 
with equal mass dijet systems. The preselection remains unchanged with respect to Ref. ||. 

A five-constraint kinematic fit is performed with energy-momentum conservation and 
equal dijet-mass constraints. In this fitting procedure, the errors on the jet energies and 



angles are parametrized as for the W mass measurement in the four-jet channel ||19|| . The 
pairing is chosen as the dijet combination giving the minimum \ 2 . 

To evaluate the mass difference between the two dijet invariant masses, momentum 
and energy conservation is imposed to rescale the energies of the four jets, fixing the jet 
velocities at their measured values. The mass difference Am between the two rescaled dijets 
is required to be smaller than 30 GeV/c 2 . 

To improve the background rejection a linear discriminant D 3 is constructed, combining 
the following five variables: 

• the production polar angle prod , i.e. the angle between the Higgs boson momentum 
direction and the beam axis; 

• the difference Am between the two rescaled dijet masses; 

• the x 2 °f the 5C kinematic fit; 

• the product of the minimum jet energy E m \ n and the minimum jet-jet angle 6? qq ; 

• the logarithm of the QCD four-jet matrix element squared A4qcd • 



The optimized LDA coefficients were determined at yfs = 206 GeV with a cocktail of 
five charged Higgs boson masses ranging between 80 and 88 GeV/c 2 , leading to: 

D 3 = -0.951 cos 2 # prod - 0.0065 Am - 0.000968 xlc 
- 0.0034 {E min x # qq ) - 0.335 log 10 (.M Q cD) , 

with Am in GeV/c 2 , E min in GeV, # qq - in radians, and .Mqcd in GeV -4 . The distribution of 
D 3 is shown in Fig. |3]a. The cut was optimized for m H ±=76, 80 and 84 GeV/c 2 . Events are 
accepted if D 3 > 1.3. For m H ±=75 GeV/c 2 and B(H + — >t + z/ t )=0, a total of 101.9 events is 
expected in the data. The efficiency does not depend on ^/s. 

After the complete selection, the comparison between data and simulation is displayed 
in Fig. |3]b for the dijet invariant mass. The numbers of events observed in the data are 
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Figure 3: (a) The distribution of the discriminating variable D%. (b) The distribution of the 
reconstructed mass of the Higgs boson candidates after the cut on the discriminating variable. The 
points are the data, the open histograms are the Standard Model backgrounds and the hatched 
histogram represents the Higgs signal expectation for m H ± = 75GeV/c 2 . The signal is arbitrarily 
normalized. 



compared in Table [I] to the expected background from Standard Model processes, dominated 
by W + W~ production. An overall 2.4 standard deviation deficit with respect to expectation 
is observed. It is correlated with the deficit observed in the measurement of the W + W~ 



hadronic cross section [121, which was ascribed to a statistical fluctuation. 



The systematic error on the number of expected signal events is estimated to be 2.5%. 
The main contributions are from limited sample statistics (1.3%), uncertainty on the cross 
section for charged Higgs production (2%) and accuracy of the simulation (0.5%). The 
systematic error on the expected background, dominated by W + W~ and qq production, 
is estimated to be 2.0%. The main contributions are from the simulated sample statistics 
(0.4% for W + W~ and 1.6% for qq), the uncertainty on the cross section (0.5% for W + W~ 
and 5% for qq), and the adequacy of the simulation (1.4% for W + W~ and 2.1% for qq). 



4 Results 



No evidence for a signal is observed in the data. The results of the three selections have 
been combined to set a 95% C.L. lower limit on the mass of charged Higgs bosons. 

Full background subtraction has been performed in setting the limit with the likelihood 
ratio test statistic [21]. Systematic uncertainties are taken into account according to 



Ref. [22]. To improve the sensitivity of the analysis, the charged Higgs boson mass has 
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been used as a discriminating variable for the esse and cst~u t channels. In the previous 
publications Jl|, ||], j3|, only event counting was used in the t + v t t~v t channel. In this 
analysis, the discriminant variable Dq has been introduced in the limit setting procedure. 

The result of the combination of the three analyses is shown in Fig. [|. Charged Higgs 
bosons with mass lower than 79.3GeV/c 2 are excluded at the 95% CL. independently 
of B(H + — >r + z/ r ). The corresponding expected exclusion is 77.1GeV/c 2 . For the values 
B(H+ ^t + v t ) = and 1, 95% C.L. lower limits on m H ± are set at 80.4GeV/c 2 (with 
78.2GeV/c 2 expected) and 87.8GeV/c 2 (with 89.2GeV/c 2 expected) respectively. 




Figure 4: Limit at 95% C.L. on the charged Higgs boson mass as a function of B(H + ^r + ^ r ). 
The expected (dash-dotted) and observed (solid) exclusion curves are shown for the combination 
of the three analyses, using the full 189-209 GeV data set. 

Upper limits can also be derived on the H + H~ cross section at y/s = 200 GeV, as a 
function of the Higgs boson mass, for B(H + — >t + v t )=0, 50 and 100%. To combine the 
data at different centre-of-mass energies, the limit on the cross section was extrapolated to 
200 GeV with the expected y/s dependence for the production of a charged scalar particle 
pair. The result is shown in Fig. || as a function of mn± • 
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Figure 5: Upper limits at 95% C.L. on the H + H production cross section at ^Js = 200 GeV 
for B(H+^r+z/ r ) = l (dashed line), B(H+ ^r+z/ r )=0 (dotted line) and B(H+^r+z/ r )=0.5 
(dashed-dotted line). The charged Higgs boson production cross section in the two-Higgs- 
doublet model is shown as a solid curve. 



5 Conclusions 



Pair-produced charged Higgs bosons have been searched for in the three final states 
t + v t t~v t , cst~P t and esse, with 629 pb -1 of data collected at centre-of-mass energies from 
189 to 209 GeV. No evidence for Higgs boson production was found and lower limits were 
set on m H ± as a function of B(H + ^t + u t ), within the framework of two-Higgs-doublet 
models. Assuming B(H + ^t + z/ t )+B(H + ^cs) = 1, charged Higgs bosons with mass below 
79.3GeV/c 2 are excluded at 95% C.L., independent of B(H+ ^r+z/ r ). 
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